Ground state entanglement in quantum spin chains 
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A microscopic calculation of ground state entanglement for the XY and Heisenberg models shows 
the emergence of universal scahng behavior at quantum phase transitions. Entanglement is thus 
controlled by conformal symmetry. Away from the critical point, entanglement gets saturated by 
a mass scale. Results borrowed from conformal field theory imply irreversibility of entanglement 
loss along renormalization group trajectories. Entanglement does not saturate in higher dimensions 
which appears to limit the success of the density matrix renormalization group technique. A pos- 
sible connection between majorization and renormalization group irreversibility emerges from our 
numerical analysis. 
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I. INTRODUCTION 

At zero temperature, the properties of a quantum 
many-body system are dietated by the structure of its 
ground state. The degree of complexity of this structure 
varies from system to system. It ranges from exception- 
ally simple cases -e.g. when an intense magnetic field 
aligns all the spins of a ferromagnet along its direction, 
producing a product or unentangled state- to more intri- 
cate situations where entanglement pervades the ground 
state of the system. Thus, entanglement appears nat- 
urally in low temperature quantum many-body physics, 
and it is at the core of relevant quantum phenomena, 
such as superconductivity ^ L quantum Hall effect Q and 
quantum phase transitions 31 ■ 

There are several good reasons to study entanglement 
in quantum many-body systems. On the one hand, over 
the last decade entanglement has been realized to be 
a crucial resource to process and send information in 
novel ways Q. It is, for instance, the key ingredient in 
quantum information tasks such as quantum teleporta- 
tion and superdense coding, and it also appears in most 
proposed algorithms for quantum computation This 
has triggered substantial experimental efforts toproduce 
entanglement in engineered quantum systems Con- 
sequently, one may want to investigate and character- 
ize entanglement in those systems where it appears in a 
natural way, with a view either to extract it to process 
quantum information or else to gain insight into physical 
mechanisms that can be used to entangle a large number 
of quantum systems. 

But one can also motivate these studies without re- 
ferring to potential applications of entanglement as a re- 
source in quantum information processing. The ground 
state of a typical quantum many-body system consists 
of a superposition of a huge number of product states. 
Understanding this structure is equivalent to establish- 
ing how subsystems are interrelated, which in turn is 
what determines many of the relevant properties of the 
system. In this sense, the study of multipartite entan- 



glement offers an attractive theoretical framework from 
which one may be able to go beyond customary ap- 
proaches to the physics of quantum collective phenom- 
ena 0- Most promisingly, a theory of entanglement in 
quantum many-body systems may also lead to the de- 
velopment of new numerical techniques. In particular, 
recent results @ show how to efficiently simulate certain 
quantum systems through a suitable parametrization of 
quantum superpositions. 

In this paper we present a quantitative analysis of en- 
tanglement in several one-dimensional spin models, ex- 
panding and complementing the results of Ref. . The 
models we discuss fulfill a convenient combination of re- 
quirements: they are solvable -the ground state can be 
computed by using well-known analytical and numerical 
techniques- and at the same time they successfully de- 
scribe a rich spectrum of physical phenomena, which in- 
clude ordered and disordered magnetic phases connected 
by a quantum phase transition [Sj. The paper has been 
divided into five more sections. A brief summary of them 
follows. 

The study of entanglement in a system with many par- 
ticles can be approached in several complementary ways 

11 E 111 111 El 111 III 111 III Im JJ« Ja] in sec- 
tion^ we briefly review some of them and explain the 
specific aspect on which we focus here. We consider a 
quantum spin chain in its ground state. Our aim is to 
determine the degree of entanglement between a block of 
spins and the rest of the chain, as measured by the von 
Neumann entropy of the block, and to investigate how 
this entanglement grows with the size of the block. 

Sections IIIII and HVl are devoted, respectively, to com- 
puting the entropy of a spin block for the XY model and 
for the XXZ model. The calculation is divided into two 
parts. First, we need to construct the ground state of the 
spin chain, which for general chains is a highly non-trivial 
problem. Fortunately, the XY model can be treated an- 
alytically in the limit of an infinite chain |2l |2l |23| . 
Similarly, for the XXZ model there are known techniques 
[2^ ') to easily cope with chains consisting of up to twenty 
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spins. Then, from the ground state of each model we ex- 
tract the entropy of a spin block. For the XY model 
this is shown to build down to diagonalizing a matrix 
whose dimensions grow only quadratically with the size 
of the block. In this way we compute the entanglement 
for blocks of up to several hundreds of spins. In the 
XXZ model, instead, we only consider blocks of up to 
ten qubits, but the results can already be convincingly 
interpreted as an independent confirmation of the con- 
clusions drawn from the XY model. 

In section|3wc discuss the findings of the previous two 
sections, that can be summarized as follows: 

• Off a critical point, the entanglement of a block 
of spins with the chain -a function that turns out 
to grow monotonically with the number of spins 
in the block- achieves a finite saturation value for 
sufficiently large blocks. 

• At a quantum phase transition^ instead, the entropy 
of a block of spins grows unboundedly. More specif- 
ically, the entropy for a critical chain grows log- 
arithmically in the size of the spin block, with a 
multiplicative coefficient that depends only on the 
universality class of the phase transition. That is, 
at the critical point entanglement obeys a universal 
scaling law. 

Interestingly, the behavior of the entanglement in a 
critical spin chain matches well-known results in confor- 
mal field theory, where the geometric entropy -analogous 
to the spin-block entropy, but defined in the continuum- 
has been computed for 1-1-1 dimensional theories p5ll2^ 
l28j . The geometric entropy grows logarithmically 
with the size of the interval under consideration and with 
a multiplicative constant given by the eentral charge of 
the theory. As described in section FVll a consistent pic- 
ture arises. At a critical point the large-scale behavior of 
a spin chain is universal, with the quantum phase tran- 
sition belonging to a given universality class. All long- 
range properties of the chain are then described by the 
conformal field theory associated with that universality 
class. In particular the entanglement between a large 
block of spins and the rest of the chain follows the same 
law as the geometric entropy. 

The above connection between entanglement and the 
geometric entropy of conformal theories has several im- 
plications. Previous calculations of the entropy in higher 
dimensions indicate that in 2- and 3-dimensional spin lat- 
tices the entropy of a spin block grows as the size of the 
boundary of the block — the same law holding both for 
critical and non-critical lattices. The lack of saturation 
of the entropy as a function of the size of the block ex- 
plains the failure of the DMRG technique (2^. Thus, we 
interpret in terms of entanglement why this numerical 
technique — so successful for non-critical spin chains — 
deteriorates in critical one-dimensional lattices and fails 
to work properly both for critical and non-critical chains 
and in 2- and 3-dimcnsional lattices 1301 . 



The scaling law obeyed by the entanglement of a criti- 
cal ground state implies that the greater a spin block is, 
the more disordered or mixed its density matrix. Thus, 
the entropy indicates an ordering of the reduced density 
matrices, according to how mixed they are. This order- 
ing can be further refined and shown to actually emerge 
from majorization relations between the reduced density 
matrices. 

Finallyj we translate the results related to the c- 

theorem jSlJ to quantum information. Entanglement is 
argued to decrease along renormalization group trajecto- 
ries. A number of numerical and analytical results are 
consistent with the idea that irreversibility of renormal- 
ization group flows may be rooted on a majorization or- 
dering of the vacuum density matrices along the flow. 

II. ENTANGLEMENT MEASURES IN A 
QUANTUM SPIN CHAIN 

A major difflculty in studying the entanglement in a 
many-body quantum system comes from the fact that 
the number of degrees of freedom involved grows expo- 
nentially with the number of interacting subsystems. In 
particular, the task of computing explicitly the ground 
state of a chain consisting of a large number of spins, to 
then analyze its entanglement properties, turns out to be 
very difflcult, if not insurmountable. 

Fortunately, for some speciflc spin models the ground 
state has been previously computed. One can then try 
to characterize entanglement in these models. This is 
again a rather ambitious enterprise. On the one hand, 
it entails serious computational difficulties, since the cor- 
responding ground states, when expressed in a local ba- 
sis, still involve exponentially many coefficients. On the 
other hand such characterization is also challenging from 
a conceptual viewpoint. The study of entanglement of a 
large number of particles is a relatively unexplored sub- 
ject and it has not yet even been established what aspects 
of a ground state a sensible characterization should con- 
sider. Therefore, an important part of the problem is to 
actually identify which quantities may be of interest. 

This section is devoted to describe and motivate our 
particular approach, which attempts to characterize the 
ground state of the spin chain through the spectral prop- 
erties of the reduced density matrix for a block of spins, 
and in particular through its entropy. We start by pre- 
senting some generalities and reviewing previous work. 

A. Overview of previous work 

A state l^*) S of N spins is entangled if it cannot 

be written as the tensor product of single-spin states, 

|*)^|7/;i)®|7^2)®---<»|^jv). (2.1) 

Product states depend on 0{N) parameters and are 
therefore just a subset of zero measure in the set of states 
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of N spins. A generic entangled state, when expressed 
in a local basis, depends on 0{2^) parameters. Char- 
acterizing entanglement is about identifying a reduced 
subset of parameters that are particularly relevant from 
a physical or computational point of view. 

1. Entanglement under local manipulation 

In recent years a quantitative theory of bipartite en- 
tanglement has been developed. As proposed in the pio- 
neering work by Bennett, Bernstein, Popescu and Schu- 
macher |32, this theory is based on the possibility of 
converting one entangled state |^) into another entan- 
gled state 1 5"') by applying local operations on each of 
the subsystems and communicating classically, a set of 
transformations denoted as LOCC (see for extensive 
reviews). The basic idea is that if the state j^") can be 
converted into the state l^*') by LOCC, 

I*) — I*') (2.2) 

then l^f) cannot be less entangled than j^"'), since LOCC 
can only introduce classical correlations between the sub- 
systems. Local convertibility can in this way be used to 
compare the amount of entanglement in different states. 

Following these ideas two remarkably simple character- 
izations of pure-state entanglement are possible for sys- 
tems with N = 2 subsystems: 

(z) Bennett et al. |33| showed that, in an asymptotic 
sense, any entangled state l^')^!^ of two particles A and 
B is equivalent (that is, reversibly convertible by LOCC) 
to some fraction E{'i') of an EPR state, 

^(|0)^®|0)s4-|1)^«)|1)b). (2.3) 

Here the entropy of entanglement El"^) corresponds to 
the von Neumann entropy of the reduced density matrix 
PA = trB|5')(\E'| for any one of the systems, 

E{^) = -tr{pAlog^PA). (2.4) 

(m) Nielsen [s^l subsequently showed that determinis- 
tic conversions of a single copy of \^) into \^') by LOCC 
are ruled by the majorization relation (to be introduced 
later in Eq. ()2.17(l ). More general LOCC transforma- 
tions are similarly ruled by a finite set of entanglement 
monotones [ssj . 

The optimal local manipulation of a bipartite system 
in a pure state is presently well understood. Bipartite 
pure-state entanglement can be characterized by a sin- 
gle measure E{'^) in the asymptotic regime and by a 
small set of entanglement monotones in the single-copy 
case. However, none of these results has been success- 
fully generalized to systems with N > 2 subsystems^. In 



^ The lack of complete generalizations of bipartite results to > 2 



spite of the remarkable success achieved for bipartite sys- 
tems, LOCC transformations do not seem to be a good 
a guidance to comprehensively characterize multipartite 
entanglement^. Nevertheless, the above results can still 
be used to characterize any bipartite aspects of the en- 
tanglement in a multipartite system, and therefore will 
be sufficing for the purposes of this paper. 



2. Entanglement in spin chains 

The study of entanglement in condensed matter sys- 
tems was initiated by Nielsen He originally ana- 
lyzed two interacting spins in the Heisenberg model with 
an external magnetic field and studied how entanglement 
depends on the temperature and the intensity of the spin- 
spin interaction and magnetic field. In his calculations, 
Nielsen used Wooters' concurrence [sH, a measure of 
mixed-state entanglement defined for two-qubit systems. 

More recently, several other authors have also stud- 
ied the concurrence in spin systems. Wooters has 
studied the maximal nearest neighbor concurrence that 
an infinite, translationally invariant spin chain can have, 
a result extended by O'Connor and Wooters ^3 to fi- 
nite spin rings. Furthermore, concurrence in the two- 
spin Heisenberg model has been reanalyzed by Arnesen, 
Bose and Vedral Gunlycke, Bose, Kendon and Ve- 
dral |14| have considered a ring of several spins with Ising 
interaction and external magnetic field, and studied how 
two-spin entanglement depended on the orientation of 



subsystems has a simple explanation, at least for single-copy con- 
versions. The most allowing scenario for local manipulation of 
entanglement in the single-copy regime is that of stochastic local 
operations SLO, where the conversion of Eq. ST^ is only re- 
quired to succeed with some non-vanishing probability. Clearly, 
if a conversion is not possible by SLO, then it is also not possi- 
ble by LOCC. But for systems with N > 3 subsystems (and with 
the exceptional case of three qubits) , two randoml y ch osen states 
I "I") and \^') are generically unconnected by SLO l36l . This can 
be understood by noticing that the total number of parameters 
accessible to local manipulation grows linearly with the number 
TV of subsystems (the most general SLO operation can be imple- 
mented by a single measurement on each subsystem), whereas 
I'l') and I*!"') depend on exponentially many parameters. 
^ There are many other possibilities to be considered instead, in- 
volving a coarse-grained look at entanglement. One could base 
a characterization of entanglement in quantum many-body sys- 
tems on the operational complexity of preparing a quantum state 
I't) (or a series of quantum states {I'fjv)} involving an increasing 
number of particles A'^) by only two-particle unitary operations. 
The entanglement of states {'^ m) B.nd l^'^y) would be compa- 
rable if, say, it only takes poly(7V) two-particle operations to 
interconvert them. Parameter counting shows that one can then 
distinguish between the class of states that can be produced with 
poly(A'^) two-particle operations and, for instance, those requir- 
ing exp(Af) operations. Alternatively, one can study the compu- 
tational cost of classically simulati ng t he state of a many-body 
quantum system and its dynamics 1371 . A third possibility is to 
consider how entanglement is affected by a change of scale in the 
system. 
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the magnetic field. Wang, Fu, Solomon [l5j have studied 
the anisotropic Heisenberg model with three spins. For a 
Heisenberg ring of N spins, Wang and Zanardi jl^ have 
expressed the nearest neighbor concurrence in terms of 
the internal energy of the ring and analyzed the violation 
of Bell inequalities, and Wang has investigated the 
concurrence in the XX model. Ref. contains some 
other related works. 

Osterlohj Amico, Falci and Fazio and Osborne and 
Nielsen have recently studied entanglement in the 
ground state of an infinite XY and Ising spin chains and 
its relation to quantum phase transitions. More specifi- 
cally, they have computed the concurrence between pairs 
of spins, for different choices of the pair. For the Ising 
model with transverse magnetic field, the concurrence at 
nearest and next to nearest neighbors has a maximum 
near the critical point. Osterloh et al. have suggestively 
noticed that there seems to be some form of universal 
scaling in the derivative of the concurrence of nearest 
neighboring spins, and similarly for the second deriva- 
tive of the next to nearest neighbors. A notable fact is 
that the concurrence seems to disappear at third nearest 
neighbors. As pointed out by Osborne and Nielsen, this 
can be interpreted in terms of the monogamy of entan- 
glement 113 ■ In practice, it can also be understood as 
a shortcoming for using a two-qubit measure in order to 
capture the global distribution of entanglement along the 
chain. 



B. Entropy of a block of spins 

The approach we follow here has been proposed by Vi- 
dal, Latorre, Rico and Kitaev and, as previous works 
based on the concurrence, it is focussed on bipartite en- 
tanglement. But instead of analyzing the entanglement 
between two of the spins of the system, we consider a 
whole block of adjacent spins and study its entanglement 
with the rest of the chain. We are particularly interested 
in how the entanglement between the block and the chain 
depends on the size of the block. In this way, we expect to 
be able to explore the behavior of quantum correlations 
at different length scales and to capture the emergence 
of universal scaling at a quantum critical point. Later 
on in this section we will further motivate our choice by 
referring to the relationship between entanglement and 
the efficiency of numerical schemes for the simulation of 
spin chains. 

Let I'^g) denote the ground state of a chain of N spins 



and let p^, 



PL = tlN-L\'^g){'^g 



(2.5) 



be the reduced density matrix for L contiguous spins. 
In the models we shall discuss, the ground state jvl/g) is 
translationally invariant, so that pL does not depend on 
the position of the block of spins but only on its size L. 
Because the chain is in a pure state, all the information 
about the entanglement between the block of spins and 



the rest of the chain is contained in the eigenvalues of 
Pl. If Pi is a pure state itself, then the block is unentan- 
gled from the chain. Instead, if pL has many non-zero 
eigenvalues, this roughly indicates a lot of entanglement 
between the block and the chain. As a matter of fact, the 
whole spectrum Sp(pl) is of interest, as we shall discuss 
shortly. For concretcness, however, we will mainly use 
a single function of the spectrum of pi, namely the von 
Neumann entropy'^ 



Sl = -tr (pi log2 pl) 



(2.6) 



as a measure of entanglement. This choice corresponds 
to the entropy of entanglement [s^ (recall Eq. H2.4|) ') 
between the block and the rest of the chain. 

L spins 

I ^ 




/ \ \ / \ 




Figure 1: The entropy Sl corresponds to the von Neumann 
entropy of the reduced density matrix pL for a block of L 
adjacent spins, and measures the entanglement between the 
block and the rest of the chain. State pL is obtained from the 
ground state |$g) of the N spin chain by tracing out all TV — i 
spins outside the block. 



1. Properties of tfie entropy of a block of spins 

Let us discuss some general properties of Sl as a func- 
tion of L. Sl is positive by construction. 



Sl>0, 



i = 0,1,' 



(2.7) 



where for convenience we define 5*0 = 0. Because the 
chain is in a pure state the spectrum Sp(pl) of the 
reduced density matrix for a block of spins and the spec- 
trum for the rest of the chain are the same. In particular, 
the two parts will also have the same entropy. Recalling 
that the ground state is translationally invariant, 
we have 



Sl^Sn-l, L = 0,1,---,N. 
In addition, Sl is a concave function |40l |. 

r, , Sl-m + Sl + m 

JL > , 



(2., 



(2.9) 



Using the logarithm to base 2, the entropy is measured in units 
of information or bits. 
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where L = 0, - ■ ■ ,N, and M = 0, • • • ,min {N-L,L}. 
This can be proved with the help of the strong subaddi- 
tivity of the von Neumann entropy 0, ^| , 



S{ABC) + S{B) < S{AB) + S{BC), 



(2.10) 



where A, B and C are three subsystems and, say, S{AB) 
denotes the entropy of pab, the joint state of systems A 
and B. Let A, B and C correspond to three adjacent 
blocks of our translational invariant spin chain, with M , 
L — M and M spins respectively. Then we have 



S{ABC) 
S{AB) = S{BC) 
S{B) 

and Eq. (|2.10|l reads 



Sl+m 
Sl 

Sl~m, 



S 



L+M 



5, 



L-M 



<2Sr 



(2.11) 
(2.12) 
(2.13) 



(2.14) 



from where Eq. H2.9|l follows. 

Finally, we note that the above properties imply that 
Sl does not decrease as a function of L in the interval 
L G [0,N/2]. In particular, in the limit of an infinite 
chain, N —f oo, Sl becomes a non-decreasing, concave 
function for all finite values of L. 



2. Examples 

It is possible to get extra insight into the properties 
of Sl as a measure of entanglement by analyzing some 
particular cases, as illustrated in Fig. ||2Jl. We note first 
that 5l is upper bounded by 



Sl < mm{L,N-L}, 



(2.15) 



since pL is supported in a local space of dimension = 
min{2^, 2^~^}, whereas Sl vanishes for all L only for 
product (i.e. unentangled) states. 

The paradigmatic GHZ state of N spins or qubits. 



II 



(2.16) 



is often regarded as a maximally entangled state. How- 
ever, from the present perspective it is only slightly en- 
tangled. Indeed, the entropies of a block of spins are 
Sl = 1 for L = 1, • • • , A^— 1, and are therefore far below 
the upper bound H2.15|l . 

Here we will be concerned with the ground state \'^g) of 
spin chains that are invariant under discrete translations 
by any number of sites. [For finite chains, we will assume 
that the extremal spins are connected (spin rings) and 
will require invariance under circular translations.] One 
could expect that translational symmetry of implies 
a more restrictive bound for the values Sl can achieve. 
However this is not the case, since Stelmachovic et al. 
[4^ have found a translationally invariant state that sat- 
urates H2.15|l . This is in contrast with the case of states 



that are invariant under all possible permutations of the 
spins. There the dimension d^™ = L -I- 1 of the symmet- 
ric subspace leads to the upper bound S'lf™ < logj L + 1. 

Finally, at a critical point the ground state may have 
some extra symmetries. In particular, it is known that 
in the large scale limit (that is, for scales much larger 
than the distance between neighboring spins) a critical 
spin chain is conformal invariant. We will explore the 
implications of this additional symmetry in section IVII 
The ground state entropy for critical chains will turn out 
to grow as Sl = fclog2 L for some universal constant k. 




Figure 2: Bounds for the entropy Sl for some pure states in 
a system with N=26 spins. Triangles correspond to the lin- 
ear upper bound 12.151 1. and applies to translationally invari- 
ant states (and, more generally, to arbitrary A'^-qubit states). 
Stars are the logarithmic upper bound for a symmetric state 
under permutations. The diamonds are the values of the en- 
tropy for a GHZ state. 



3. Majorization and von Neumann entropy 

As mentioned above, the entanglement between a block 
of spins and the rest of the chain is a function of the spec- 
trum Sp(/9l) of the reduced density matrix of the block. 
Our ultimate aim is to characterize how this entangle- 
ment depends on the number L of spins in the block. 
A main motivation for this is that in this way we hope 
to capture the emergence of universal scaling for entan- 
glement at a quantum phase transition. Therefore we 
would like to be able to compare the spectrum Sp(pi,) 
for different values of L. 

The entropy Sl can be used for this purpose, for it 
establishes an order in the set of probability distributions 
-equivalently, in the set of spectra of density matrices. 
For instance, we have mentioned above that the entropy 
Sl is non-decreasing in the interval L e [0, A^/2]. We can 
now use this result to say that, according to the entropy, 
the entanglement of a block of spins and the rest of the 
chain monotonically increases with the size of the block 
(for blocks smaller than half of the chain). 

Nevertheless, there are other powerful tools to com- 
pare probability distributions, and by using them one 
may obtain a finer characterization of entanglement. In 
particular, a far more tight sense of (partial) ordering 
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between probability distributions is established by the 
majorization relation |43l | , a set of inequalities that con- 
trol the conversion of bipartite entanglement by LOCC 
in the single-copy scenario [3 HE Q • 

Let us briefly recall that a given probability distri- 
bution X = {xi} (where xi > ■ ■ ■ > Xn) is majorized 
by another probability distribution y = {yi} (where 
2/1 > • • • > J/n), denoted x ^ y, when the following scries 
of inequalities arc simultaneously fulfilled: 

xi < vi 
X1 + X2 < yi + y2 

xi+x2-\ \-Xn = yi+ij2-\ \-y,i- (2.17) 

The majorization relation x < y expresses the fact that y 
is more ordered than x. Given two arbitrary probability 
distributions x and y, inequalities H2.17|l are not likely to 
be simultaneously fulfilled, but when they are, most mea- 
sures of order are consistent with x < y. In particular, 
the von Neumann entropy fulfills 

p<p' ^ S{p)>S{p'), (2.18) 

where p < p' refers to majorization between the spectra 
of these two density matrices. 

In section IVTI we shall explore whether a majorization 
relation underlies the scaling behavior of the entropy Sl 
at critical points. 



C. Entanglement 
in numerical studies of a quantum spin chain 

There are many aspects of the ground state of a spin 
chain that could be taken as a guide to characterize its 
entanglement. Our choice can be motivated by the role 
the reduced density matrix of a block of spins plays 
in some numerical schemes. We finish this section by 
explaining how the spectrum Sp(pl) of pL determines 
the efficiency of White's density matrix renormalization 
group (DMRG) method and of a recently proposed 
simulation scheme for the simulation of quantum spin 
chains 0. 

White's DMRG method is a numerical technique 
that has brought an enormous progress in the study of 
one-dimensional quantum systems such as quantum spin 
chains. It allows to compute ground state energies and 
correlation functions with spectacular accuracy for non- 
critical spin chains. The DMRG method, however, loses 
its grip (as many other methods) near a critical point 
and fails to work for quantum spin lattices in two or 
three dimensions even away from the critical poin t [30l |. 
As recently explained by Osborne and Nielsen [43 in the 
language of quantum information, the degree of perfor- 
mance of this method is directly related to the way it 
accounts for entanglement. 



Let us consider a large spin chain and its ground state 
Y^g). The DMRG method is based on computing prop- 
erties of \^g) by constructing an approximation to the 
reduced density matrix p^ for a block of L spins, for an 
increasing value of L. This is done by retaining only the 
relevant degrees of freedom of the Hilbert space associ- 
ated to the block of spins. Such degrees of freedom are 
given by the eigenvectors of pL with greatest weights or 
eigenvalues {p\}- that we assume dccreasingly ordered, 
p'l>pT (* = 1,--- ,2^-1). 

Notice that the spectrum Sp(pl) typically contains as 
many as 2^ relevant eigenvalues, in which case the com- 
putational cost of the DMRG method explodes as the 
block size L grows. However, not all eigenvalues have 
the same weight, and a good approximation to p^ may 
be significantly cheaper to achieve than the exact reduced 
density matrix. Let denote the number of eigenvalues 
such that 

XI 

Y.pl>l~e, 0<e«l. (2.19) 

1=1 

That is, is an effective rank of pL, resulting from 
ignoring all the smallest eigenvalues that sum up less than 
e. Then, if we are willing to accept a degree of accuracy 
e, the DMRG method need only retain Xi eigenvectors 
of PL- 

The efficiency of the DMRG depends on how small Xl 
is. In turn, X'l depends on how fast the eigenvalues 
decay with i or, relatedly, on the entanglement between 
the block of spins and the rest of the chain. A very spread 
spectrum, roughly equivalent to a lot of entanglement, 
translates into a large and a large computational cost. 
If, instead, the eigenvalues {p}^} decay very fast with i, 
implying that there is not much entanglement, then X)^ 
is small and so is the computational cost of the DMRG 
method. 

Therefore, by studying the spectrum of pL, and in par- 
ticular the effective rank Xi, we may be able to assess 
how well the DMRG will perform for given values of the 
parameters (external magnetic field, spin-spin interac- 
tion) defining a particular spin model. 

On the other hand the effective rank Xi appears also 
as a decisive parameter in a recently proposed numeri- 
cal scheme for the classical simulation of quantum spin 
chains 0- In this scheme the cost of the simulation is 
linear in the number N of spins in the chain and grows 
as a small polynomial in X*^, 

X'=maxXl, (2.20) 

that is, polynomial in the maximal effective rank achieved 
for blocks of adjacent spins. 

Summarizing, the spectrum of pi, through the effec- 
tive rank Xi , is of direct interest for the numerical study 
of spin chains. The entropy Sl of pl is related to the 
effective rank of pl- Indeed, we have 

Xl""°>2^^. (2.21) 
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In addition, numerical evidence in spin chains indicates 
that 2'^^ also gives a rough estimate of for small e > 0. 

In section Owe shall discuss the results we have ob- 
tained for the entropy Sl, both for critical and non- 
critical spin chains, and analogous results for spin lat- 
tices. We will conclude that the degree of performance 
of the DMRG method for spin systems depends on how 
the entanglement between a block of spins and the rest 
of the system scales with the size of the block. 

Note added: after completing the present work we have 
become aware of a number of contributions by Peschel et 
al ^3 that study the spectrum of the reduced density 
matrix also with a view to assess the performance of 
the DMRG. 



III. XY MODEL 



In this section we study the entanglement of an infinite 
XY spin chain. We start by reviewing the main features 
of the XY model and by identifying some of the critical 
regions in the space of parameters that define the model. 
Then, we proceed to compute the ground state |^'g) of 
the system, from which we obtain the reduced density 
matrix for L contiguous spins. The knowledge of the 
eigenvalues of pi allows us to compute its entropy Sl 
and; therefore, have a quantification of entanglement in 
spin chains. Further information contained in the eigen- 
values of the density matrix will be explored in section 

Ell 

The calculations of the spectrum and ground state of 
the XY model that appear in this section and in the 
appendices review previous work in spin chains. Lieb, 
Schultz and Mattis solved exactly the XY model 
without magnetic field; Katsura computed the spec- 
trum of the XY model with magnetic field; Barouch and 
McCoy obtained the correlation function for this 
model. Finally, the entropy Sl was computed by Vi- 
dal, Latorre. Rico and Kitaev |^. Here we shall present 
an expanded version of this computation. 



field, as given by the Hamiltonian^ 

(3.1) 

Here I labels the N spins, (t[' (p = x, y, z) are the Pauli 
matrices, 

(3.2) 

acting on spin /, with 

K-'^m] = 2i5,™ ^ W^r> (3-3) 

r—x^y^z 

whereas parameter A is the intensity of the magnetic field, 
applied in the z direction, and parameter 7 determines 
the degree of anisotropy of spin-spin interaction, which 
is restricted to the xy plane in spin space. 

The XY model encompasses two other well-known spin 
models. If the interaction is restricted to the x direction 
in spin space, that is 7 = 1, then Hxy turns into the 
Ising Hamiltonian with transverse magnetic field, 

i?isi„g = -^5](araiVi+Aaf). (3.4) 

I 

If, instead, we consider the interaction to be isotropic in 
the xy plane, 7 = 0, then we recover the XX Hamiltonian 
with transverse magnetic field, 

Hxx = -\Y^ (^K^f+i + <J + Aaf ) . (3.5) 

We note that these Hamiltonians are used to model 
the physics of one dimensional arrays of spins, but also 
to describe other quantum phenomena. For instance, the 
XX Hamiltonian corresponds to a particular limit of the 
boson Hubbard model 




where a are bosonic annihilation operators, 

[aua\^] = 5i„,, (3.6) 

and ni = cL^ai are number operators. The boson Hubbard 
model (see chapters 10 and 11 of Q) consists of spinless 



A. The XY Hamiltonian 



The XY model consists of a chain of N spins with 
nearest neighbor interactions and an external magnetic 



* Notice that the sign of the interaction can be changed by ap- 
plying a 180 degree rotation along the z axis (in spin space) 
to every second spin. Since this is a local transformation, the 
ground state of the original and transformed XY Hamiltonians 
are related by local unitary operations. Therefore both ground 
states are equivalent as far as entanglement properties are con- 
cerned. In this sense entanglement depends on fewer details than 
other properties of the chain, such as the magnetization. 
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bosons on N sites, representing, say, Cooper pairs of elec- 
trons undergoing Josephson tunneling between supercon- 
ducting islands or helium atoms moving on a substrate. 
The first term in Hb, proportional to w, allows hopping 
of bosons from site to site. The second term determines 
the total number of bosons in the model, with fi the 
chemical potential. The last term, with L/ > 0, is a re- 
pulsive on-site interaction between bosons. Now, in the 
limit of large [/, no more than 1 boson will be present 
at each site. Thus, each of the TV sites has an effective 
two-dimensional local space, and the identification 



ai + a[, 
-i{ai - a]), 
1 - 2a|ai, 



(3.7) 
(3.8) 
(3.9) 



takes Hb into 



= -E (f [^r^r+i + ^r^f+i] - ^^r) , (3.10) 

which is the Hamiltonian of the XX spin chain model 
with transverse magnetic field of Eq (|3.5|) . 



B. Spectrum of Hxy and critical properties 

In appendix 1X1 we show that the spectrum of Hamilto- 
nian Hxy in Eq. (|3.1|) is given, in the limit of large N, 
by 



= -\/ (A - cos 0)2 + 72 sin^ 1 



(3.11) 



where </> € [— tt, tt] is a label in momentum space. This 
result was obtained by Katsura |22J . 

We can use the explicit expression ()3.11|l of to dis- 
cuss the appearance of critical behavior in the XY model 
as a function of parameters (7, A). 



The correlation length ^ characterizes the exponential 
decay of correlations in the spin chain [4^ l48l| , 

- K) «l) - cxp(-L/C), (3.12) 

whereas the low energy dispersion A is given by 

A = A^=o- (3.13) 

The critical scaling of these two quantities is character- 
ized in terms of |A — Ac| (deviation from the critical mag- 
netic field Ac) and critical exponents v and s through 



A^ |A- Ae^- 



(3.14) 



In addition, the dynamical critical behavior is given by 
the energy dispersion, 



A^^o~r(l + (0C)""), 



(3.15) 




Figure 3: Energy of the system for different values of the 
parameters A and 7 as a function of 0. The thick plot corre- 
sponds to the XX model without magnetic field, the dashed 
one to a system with A = 7 = 0.5, the flat one is the Ising 
limit without magnetic field and the dot-dashed one and the 
thick dashed plot correspond to the isotropic and Ising model 
with A = 1, respectively. 



where z is the dynamical exponent. An analysis of the 
scaling phenomena gives the following useful relation be- 
tween the critical exponents: z = sjv. 

For A = 1 and any value of the anisotropy 7 e [0, 1], 
the spectrum A^ in Eq. H3.11() has no mass gap, A = 0, 
and the spin chain is critical, so that Ac = 1. As far as 
criticality is concerned, we need to distinguish two cases 
depending on the anisotropy 7. 

(i) For 7 g (0, 1], Eq. I|3.11|l imphes that the behavior 
for the energy dispersions are 



A 



|A-1|, 



1 + ^ 



(3.16) 



At the critical point Ac = 1, we find the critical expo- 
nents z = 1 and s = 1, whereas the divergence of the 
correlation length in this interval is 



1 



lA-ll 



(3.17) 



with a critical exponent v = \. For later reference, we 
state that in this case the quantum spin chain belongs to 
the same universality class as the classical Ising model in 
two dimensions (or quantum Ising model in one dimen- 
sion). The critical behavior in this class is described by 
the conformal field theory of a free masslcss fermion in 
1+1 dimensions, with central charge equal to 1/2. 
(m) For the case 7 = 0, we have 



Arf, = I A — cos ( 



(3.18) 



and the long wave dispersion and the energy dispersion 
are 



A = A^=o = |A-1| 

A^^o- IA-I + 0V2I 
A- 1 



1 



(3.19) 
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Therefore the critical point Ac = 1 leads to the critical 
exponent z = 2 and s = 1 and the divergence of the 
correlation length in this interval is 



H 



XY 



Jordan- 
Wigner 



A 



Fourier ■ 
Bogoliubov 



B 



Hamiltonian of 
the XY model 



(3.20) 



with a critical exponent = ^- Notice, however, that in 
this case the spectrum is gapless for any A S [0, 1], 
since continuously vanishes for = arccos(A). This 
implies that the spin chain is actually critical for any 
value A e [0, 1] of the magnetic field. Again for later 
reference, we mention that in this case the quantum spin 
chain belongs to the universality class described by a free 
massless boson in 1+1 dimensions. This conformal the- 
ory has central charge equal to 1. 



critical 
Ising 




Figure 4: Some critical regions in the parameter space (7, A) 
for the XY model. The Ising model, 7 = 1, has a critical 
point at A = 1. The XX model, 7 = 0, is critical in the 
interval A € [0,1]. The whole line A = 1 is also critical. A 
complete analysis of the critical regions in this model was 
done by Barouch and McCoy in 



local 
diagonalization 



^ g Correlation matrix for 

X the whole chain 



Correlation matrix for 
a block of L spins 



Entropy of a block 
v3j^ of L spins 

Figure 5: This road map describes the steps followed in order 
to obtain the entropy Sl of L contiguous spins from an infinite 
XY chain. We diagonalize the Hamiltonian Hxy by rewriting 
it first in terms of Majorana operators a and then in terms of 
Majorana operators b. The ground state \'^g} is characterized 
by a correlation matrix F^ for operators b, for operators a. 
Correlation matrix F^ describes the reduced density matrix 
pL for a block of L spins. Sl is finally obtained from r£, the 
block-diagonal form of F^. 



C. The ground state 

We now turn to determine the ground state j^'g) of the 
XY model with open boundary conditions. 



Hxy — — 



2 I i+i 2 



(3.21) 



Summarizing, by analyzing the spectrum of Hxy one 
finds two distinct critical regions in the parameter space 
(7, A), namely the line Ac = 1 and the segment (7, A) = 
(1, [0, 1])^. The critical XX model corresponds to an un- 
stable fixed point with respect to the anisotropy 7. If 
we depart from (7, A) = (0,1) by a small perturbation 
7 7^ 0, the critical behavior of the spin chain turns from 
the universality class of the XX model into that of the 
Ising model. 



^ Barouch and McCoy 1231 have shown that also in the line defined 
by ■y^ + \^ = 1 two-spin correlators decay as a power of the 
distance between spins, the signature of criticality. 



in the limiting case of an infinite chain, N 00. 

Through a Jordan- Wigner transformation, Hamilto- 
nian Hxy can be cast into a quadratic form of fermionic 
operators, which in turn can be diagonalized by means 
of two additional canonical transformations, namely a 
Fourier transformation and a Bogoliubov transformation 
(see Appendix lAlfor details). Next we will determine the 
ground state \'9g) through a more convenient — although 
essentially equivalent — procedure that uses Majorana 
operators instead of fermionic operators. 

The present calculation was sketched in ^ and uses 
the formalism described in Originally, the ground 

state of the XY model was determined by Licb, Schultz 
and Mattis |23] in the case of no magnetic field, and by 
Barouch and McCoy in the case of magnetic field. 
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1. Majorana operators 

For each site I of the iV-spin chain, we consider two 
Majorana operators, 0,21-1 smd 0,21, defined by 



(N-l)/2 



2 



(N-l)/2 



spin 
sites 



\m<; / 



fn<)^f- (3-22) 

\m<l ) 



j-N, 
-N+1 , 



N-2 J labels for 
] Majorana 
^ operators 



Operators are Hermitian and obey anti-commutation 
relations, 



(3.23) 



The change of variables of Eq. (|3.22(1 . paraUel to the 
Jordan- Wigncr transformation described in Appendix [XI 
implies 



Figure 6: Through transformation 13.2211 . we can associate 
two Majorana operators, a2i-i and d2i, to site I of the spin 
chain. Notice, however, the non-local character of such trans- 
formation: 0,21-1 and a,2i are a product of Pauli matrices from 
sites -(A^ - l)/2 to L 



be a new set of Majorana operators, 



h\ = bp, 



{bp,bg}^2Spg. (3.33) 



a2ia2i+i = afaiaf^i=iafaf^j^, 
0,21-10,21+2 = a-fa-fcrf^i = -iafaf_^j^, 
0,21-10,21 = o-fo-f = icy I, 
so that Hamiltonian Hxy becomes 

N-l 

1 + 7 _ _ 1 — 7 



(3.24) 
(3.25) 
(3.26) 



The canonical transformation induced by W is parallel to 
the Fourier and Bogoliubov transformations for fermionic 
operators that appear in Appendix ^ where also an ex- 
plicit expression for A^ is displayed. In terms of operators 
6, Hxy reads 



H 



XY 



E 



i=~ 



0,21021 + 1 



021-1021+2 



H 



. N-l 



XY 



p^q=-N 



(3.34) 



or, equivalently, 
H 



1O21 , 



N-l 



(3.27) 



N-l 
2 

- y 

4^ 



Ak{b2k-ib2k - &2fe^2fc-i)-(3.35) 



XY 



(3.28) 



■m.7i— — N 

where A is a real, skew-symmetric matrix given by 



A: 



Aa Ai 
-A\ Ao Ai 



-A\ Ao Ai 



(3.29) 



2,. Correlation matrix 

The diagonalization of Hxy essentially concludes with 
the determination of the explicit form of matrix W , as 
presented in Appendix ^ However, in order to ana- 
lyze the resulting ground state |\E'g), it is convenient 
to momentarily switch to the more familiar language 
of fermionic operators. We define a set of A'' spinless 
fermionic operators 6, 



6fc 



b2k-i + ib: 



'2fc 



(3.36) 



and 



2A 

-2A 



, Al 



-(1-7) 
1+7 



(3.30) 



-(A^- l)/2 < fc < (Af - l)/2, obeying the anticommuta- 
tion relations 



Let W £ SO{2N) be a special orthogonal matrix that 
brings A into its block diagonal form B ~ WAW'^ , 



{bl,bp}^Skp, {bk,bp} = 0, 



(3.37) 



in terms of which Hamiltonian Hxy becomes, up to an 
irrelevant constant. 



1 

-1 



(3.31) 



Hxy = Afcfo^Sfc. 



(3.38) 



and let 



fe=- 



N-l 



bp= J2 ^prno-m, -N+l<p<N, (3.32) 



m=-N 



The ground state of Hxy is annhilated by all b, 

bk\^g)=0, (3.39) 
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so that {'^ g\b\bk\'i' g) — that is, the expectation value of 
a positive operator — vanishes and {'i/g\HxY\'^g) = 
corresponds to the smahest eigenvalue of Hxy- Then, 
since blbk + bkb\. 



I and &[&fc|*g) 



0. we also have 



bkbl\^ 



(3.40) 



Let {M) denote the expectation value (4'g|A/|4'g) for 
operator M . We readily have 



ih) = 0, 
{bkbp) = 0, 
{bkbl) = Skp. 



(3.41) 
(3.42) 
(3.43) 



More generally, Wick's theorem establishes that any non- 
vanishing expectation value corresponding to a product 
of operators b and b^ can be expressed in terms of {bkbj,) 

and (bkbp) and their complex conjugates. For instance, 
we have 



(KKblbl) = {bM{b'kAJ-M,)M 



(3.44) 



This means that is a gaussian state, completely char- 
acterized by the expectation values of the first and second 
moments, Eqs. (|3.41ll - (|3.43|l . 

We can now return to the Majorana operators b. An 
equivalent characterization of is given in terms of 
the correlation matrix {bpbq) = 6pq + *r^g, where 



k=- 



1 

-1 



(3.45) 



As direct substitution shows, amounts for both 
the expectation values (bkbj,) and (bkbp) simultaneously, 
which is the ultimate reason to conduct the present 
derivation in terms of Majorana operators. 

Finally, we use F^ to obtain the correlation matrix 
{hmO-n) = ^m,n + i^mn '^^ ^hc Original Majorana operators 
a, where F'^ = W^T^W. As shown in Appendix IbI one 
obtains 



Ho Hi 
-Hi Ho 



n 



N-1 



n, 



gi 

-g-i 



-Hn-i Ho 

(3.46) 

with real coefficients gi as given, in the limit of an infinite 
chain, N —^ oo, by 



1 

ZTT 



2tt 



-il4> 



coscj) — A — i^sin(, 
I cos (j) ~ X — i"f sin ( 



(3.47) 



We conclude that Eqs. H3.46|l - H3.47|l contain a com- 
plete characterization of the ground state j'l'g) of Hxy- 



D. Entropy of a block of spins 

The entropy of the reduced density matrix p for L ad- 
jacent spins. 



Sr 



-tr(plog2 p), 



(3.48) 



can be computed from F'^, Eq. H3.46|l . as follows. 

In the limit of an infinite chain, the middle of the chain 
is fully translational invariant, in that the same pL de- 
scribes the state of any block of L contiguous spins. For 
notational convenience we choose the block to contain 
qubits / = 1, ■ ■ ■ , i. We can expand the density matrix 
PL oi the block as 



PL 



2 

Ml 



Pfj-i-'-fj-L 

.p.L—0.x,y,z 



'L ' 



where coefficients p^-^^---^^ are given by 



(3.49) 



(3.50) 



In spite of the non-local character of transformation 
(|3.22() . the density matrix p^ can be reconstructed from 
the restricted 2L x 2L correlation matrix 



where 



-*(rf). 



m, n 



I,-- - ,2i, (3.51) 



Ho Hi 
-Hi Ho 



n 



L-l 



-n 



L-l 



Indeed, the symmetry 



n 



H 



XY 



Ho 



= H 



(3.52) 



XY 



(3.53) 



implies that p^^i-..^^^ — whenever the sum of p^s equal 
to X and of /i's equal to y is odd. For instance, for i = 4, 
terms such as poxOz and PxyOy vanish. Therefore non- 
vanishing coefficients p^j...^^^ correspond to the expec- 
tation value of a product of Pauli matrices with an even 
total number of ax^s and CTy 's. Such products are mapped 
through the inverse of transformation H3.22|l into a prod- 
uct of an even number of Majorana operators dm, with 
m € [1,2L]. (See Eqs. for an example). We 

can then use Wick's theorem to express such products in 
terms of the second moments {hmdn), m,n & [\,2L\, all 
of which are contained in F;^. 

In principle, then, one could use F;^, Wick's theorem 
and the inverse of transformation H3.22|l to compute p^, 
and extract from its spectral decomposition. How- 
ever, the spectrum of pL, and its entropy Sl, can be 
computed in a more direct way from F;^. 
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Let V e SO{2L) be such that it brings into its 
block-diagonal form r£ = VTf^V'^, 



1=1 



lyi 
-VI 



(3.54) 



Matrix V defines a set of 2L Majorana operators 



2L 
n=l 



with correlation matrix (cmCn) given by 



L Jmn- 



(3.55) 



(3.56) 



The structure of r£ implies that mode C2i-i is only cor- 
related to mode C21 , a most convenient fact that we next 
exploit. 

Again for the sake of clarity, we complete the present 
reasoning using a more familiar language of fermionic 
modes. We define L spinless fermionic operators 



ci 



C2/-I + «C2/ 



{ci,C,n} = 0, {cl,C,n}=Sl„ 

By construction they fulfill 

c 1 + t'm 



0, 



(3.57) 
(3.58) 

(3.59) 



which means that the L fermionic modes are uncorre- 
lated, that is in a product state, 



PL = gi® ' 



QL- 



(3.60) 



[Notice that this tensor product structure does not cor- 
respond in general to a factorization into local Hilbert 
spaces for the L spins, but is instead a rather non-local 
structure]. The density matrix gi has eigenvalues 



and entropy 



S{qi) = H2 



1^1 



(3.61) 



(3.62) 



where H2{x) = —a; log a; — (1 — .t) log(l — x) denotes the 
binary entropy. The spectrum of results now from the 
i-fold product of the spectra of the density matrices gi, 
and the entropy of is the sum of entropies of the L 
uncorrelated modes, 



St 



L 

E 

1=1 



Ho 



l + vi 



(3.63) 



Summarizing: for arbitrary values of the anisotropy 7 
and magnetic field A, and in the thermodynamic limit 



corresponding to an infinite chain {N 00), the entropy 
Sl of the ground state of the XY model can in practice 
be obtained by (i) evaluating Eq. (|3.47() numerically for 
Z = 0,---,L— 1, (m) diagonalizing in Eq. H3.52|l . 
so as to obtain i/m, and (iii) evaluating Sl using Eq. 
(|3.63(l . Appendix El contains an analytical expression of 
the coefficients gi in Eq. H3.47|) for several particular 
cases. 



Scaling of the entropy 

We can now proceed to compute the entropy for the 
XY model with different parameters. It is first important 
to note that the actual diagonalization to be performed 
takes place in a 2L x 2L space, not in the huge 2^ x 2^ 
space associated to the vacuum density matrix. This is 
obvious in Fig. |(7J) where the computation can easily 
include hundreds of spins. 

The result obtained for the reduced density matrix of 
L spins in the isotropic XX model, 7 = 0, with no ex- 
ternal magnetic field, A = 0, perfectly fits a logarithmic 
behavior 



1 



■ log2 L - 



7 = 0, A = 0, 



(3.64) 



where a is a constant close to 7r/3. A least square fit 
gives an standard error of 3 10^® in the constant of the 
logarithmic term for this model. This result shows that 
entanglement of the vacuum state scales at this critical 
point, pervading the whole system and carefully organiz- 
ing the complicate superposition of states that will wind 
up reproducing correlators. The scaling of entanglement, 
furthermore follows some universality properties we shall 
discuss later. 



3 . 5 

3 

2 . 5 

2 

1.5 

1 

0.5 






200 



Figure 7: Entropy of the reduced density matrix for L spins in 
the isotropic XX model, 7 = 0, with different external mag- 
netic field A. The maximum entropy is reached when there 
is no applied external field. The entropy decreases while the 
magnetic field increases until A = 1 when the system reaches 
the ferromagnetic limit and the ground state is a product 
state in the spin basis. 
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It is easy to extend our computation to other critical 
and non-critical points in the parameter space for the XY 
system. Fig. (jHl shows the scaling of entanglement as we 
scan 7. Note again the logarithmic scaling of the entropy 
although its coefficient is now 1/6 instead of 1/3. The 
constant correction to the logarithmic scaling is such that 



so that 



hm [5i(7 = l)-5L(7)] = -7:log2 7 

L—>oo 



(3.65) 
(3.66) 




-1.4-1.2-1-0.8-0.^0.4-0.2 
LogjY 



Figure 8: Difference of the entropy 

limi_,oo [S'l(7 = 1) — 5'l(7)] for different values of the 
anisotropy 7. For every 7 the model is critical and the 
entropy scales as l/Glogj L + a{'y), where the L- independent 
function 0(7) is perfectly fitted by — | logj 7. 



When 7 = 1 the system is described by the Ising model. 
The entropy behavior in this limit can be viewed in the 
Fig. When the magnetic field is turned to A = 1 the 
entropy reproduces the scaling law 



rtlsuig 



\ log2 L 



7 = 1, A =1. 



(3.67) 



In this model, the least square fit of the logarithmic be- 
havior gives a standard error of 4 10~^. For the Ising 
model with no external field, the ground state that min- 
imizes the total energy is the Neel state, a macroscopic 
GHZ state, for which the entropy is always equal to one. 



IV. HEISENBERG MODEL 

A. The XXZ Hamiltonian 

The XXZ model consists of a chain of N spins with 
nearest neighbor interactions and an external magnetic 
field, as given by the Hamiltonian 



H 



XXZ 



E 

I 



1, 




200 



Figure 9: Entropy for the reduced density matrix to L spins 
of the Ising model, 7 = 1, with different values of the external 
magnetic field, A G {0, 1}. The maximum entropy is reached 
at the critical point when the applied field is one. Other values 
of the magnetic field lead to saturation of the entropy. For a 
magnetic field A = the ground state of the system is in the 
Neel state or GHZ state. 



As in the previous section, / labels the A'^ spins and 
[p, = x,y,z) are the Pauh matrices, Eqs. (|3.2() - (|3.3|l . 
Parameter A evaluates the anisotropy, in the z direction, 
of the antiferromagnetic Heisenberg interaction, whereas 
A is the strength of a magnetic field applied in the z 
direction. 

The XXZ model includes as special cases two other 
well-known spin models. The XXX model corresponds 
to a fully isotropic interaction, A = 1, 

Hxxx = J2 (^K^^r+i + '^f^f+i + f^z'^^f+i] + ^<^i 

(4.2) 

Also, when the interaction in restricted to the plane xy 
in spin space, A = 0, we recover the XX model of Eq. 

These Hamiltonians are commonly used to model the 
physics of certain spin chains, but also to describe other 
quantum systems. For instance, the XXX Hamiltonian 
without magnetic field can be obtained in a particular 
limit of the fermion Hubbard model, 



I 



where h are fermionic annihilation operators. 



(4.3) 



(4.4) 



(4.1) 



niT = hj^hiT are fermion number operator and I labels 
one of N sites of a chain while jy denotes one of two spin 
orientations, | or |. The fermion Hubbard Hamiltonian 
(see chapter 10 of |3|) was originally introduced to de- 
scribe the motion of electrons in transition metals, and 
consists of spin-1/2 fermionic particles moving along the 
A^ sites. Parameter e is the energy cost of having one 
fermion, t is the tunneling parameter and U quantifies 
the interaction between two fermions at the same site. 
Notice that the total number of fermions n = nir is a 
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constant of motion. Then, when restricted to a subspace 
with a given value of n, the first term in the Hamiltonian 
is proportional to the identity and can be omitted. Let 
us consider the case when the total number of fermions is 
N, that is, the same as the number of sites. In the limit 
U ^ t, two fermions are not energetically allowed to be 
on the same site, and we have one fermion per site. In 
this limit, and with the identification 



erf = h\{h^ - h\Jiii, 



(4.5) 
(4.6) 
(4.7) 



the Hubbard model can be recast (using perturbation 
theory) into the isotropic Heisenberg model without mag- 
netic field, 



I 



(4.8) 



where J = 4|j-. Metal-insulator transitions, supercon- 
ductive systems or magnetic properties can be explained 
with this Hamiltonian. 

Quantum phase transitions are identified with points 
of non-analyticity in the ground state energy of the spin 
chain. This non-analyticity may appear in the limit of 
a large chain or, as in the present case, may be due to 
level crossing. Notice that we can decompose Hxxz, Eq. 
()4.1f) . in terms of three commuting parts. 



Hxxz ^Hi + AH2 + XH3 



As we change A or A, an excited state may see its en- 
ergy decreased enough as to become the ground state. 
At one such point, the ground state energy E{A, A) will 
not be analytical. Thus, in this model quantum phase 
transitions occurs for finite chains. 



B. Bethe Ansatz for the Heisenberg model 

Several properties make the Heisenberg Hamiltonian 
with periodic boundary conditions completely integrable. 
Two symmetries are essential to get the model solution 
and are used in the Bethe Ansatz {BA) Rotational 
symmetry about the z-axis in spin space implies that the 
z-component of the total spin S^, 



(4.10) 



is conserved. Sorting the basis vectors according to the 
quantum number = N /2 — r, where N is the number 
of sites in the chain and r the number of spins down, 
is all that is required to block diagonalize the Hamilto- 
nian. The second symmetry is the invariance of H with 
respect to discrete translations by any number of lat- 
tice spacings. To reconstruct the whole spectrum of the 



Heisenberg model, Bethe's idea is to start with the fer- 
romagnetic state \F) , 



= ITT •••!>, 



(4.11) 
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whose spin angular momentum, is maximum, and to 
get a translationally invariant eigenstate of the Hamilto- 
nian with one unit less of spin angular momentum. The 
rest of the spectrum of the Heisenberg model is obtained 
by iterating this process. 

Here, we shall use the BA to get a numerical but ex- 
act solution of a finite spin chain. Finite size effects 
are present but can be controlled by comparing differ- 
ent sizes. Scaling of entanglement is thus approached 
asymptotically as the size grows. For a fixed value of 
5y = N /2 — r. eigenstates are of the form 

|*> = a{ni,...,nr)\ni,...,nr), (4.12) 

l<rii<...<?i^<Ar 

where rii, • • • , list the position of the r spin that are 
down, and a(ni, • • • , n^) fulfills 



a(ni, ...,nr) = ^ exp i ^ k^jUj + ^ ^ 6'pi,-pj 

•P6<Sr \ j = i i<j 

(4.13) 

Here V € Sr denotes one of the r! permutations of 
{l,...,r} and ki and Oi^j with (i,j) G {l,...,r} are the 
parameters to be determined. 

Three general conditions hold for these parameters 



cot 



'1,3 



A sin 



V{i,j}, 



2 " 

Nki = 2nX, + Y0^^J 



e {l,...,r}. 



ze{l,...,r}, (4.14) 



where the integers Xi are called Bethe quantum numbers. 
The states are completely determined by the Bethe quan- 
tum numbers and Eqs. H4.12|) - (|4.14|) . 

It is known [sJl that the set {Xi} with the lowest energy 
for each = N/2 — r satisfies: 

N 

X, = S^-l + 2i^ —-r-l + 2i ie{l,...,r}. (4.15) 
The expression for the ground state energy reads 

i=l 

(4.16) 

We have found convenient to solve the previous sys- 
tem of non-linear equations using a minimization code 
for absolute errors based on genetic algorithms, although 
many other numerical techniques could have been used. 
Precision can be controlled by setting a maximum bound 



15 



in the numerical error in the parameters 0's and fc's to 
10~^. This small error has proven enough to detect the 
scaling of entropy. The systems that we study have an 
even number of sites N, so the eigenvalue of the op- 
erator in the ferromagnetic system is the integer N/2. 



C. Entropy of a block of spins 

The fact that the BA is used to compute numerically 
the eigenvalues of the reduced density matrix of blocks 
of spins limits the size of the system that can be studied. 
The results we obtain are thus less precise than those for 
the XY model although scaling laws can be inferred with 
confidence. As a first step, we concentrate on the effect 
of the finite size system for the entropy results. More pre- 
cisely, we analyze the isotropic model without magnetic 
field in a chain of TV = {8, 10, 12, 14, 16, 18} sites. The 
results are plotted in Fig. (|l()|l . Finite size effects bend 
down the entropy when the size of the block approaches 
half of the chain. Smaller blocks are less sensitive to the 
finite size and show good scaling. The numerical results 
indicate that the entropy behavior converges to a loga- 
rithmic scaling when as the size of the system increases. 
This asymptotic behavior corresponds to 



Sl 



1 



log2 L 



(4.17) 



As we shall discuss in the next two sections, this entropy 
scaling falls into the universality class of a free boson. 
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Figure 10: Dependence of the entropy Sl on the finite size of 
the cliain in the isotropic Heisenberg model without magnetic 
field. In the plot, triangles, stars, diamonds, triangles, stars 
and diamonds correspond to 18, 16, 14, 12, 10, and 8 spins 
chains, respectively. All curbs coincide with the upper one for 
low block sizes L, but disagree as L gets closer to N/2 due to 
a finite-size saturation effect. 



1. Isotropic model in a magnetic field 

We first analyze the XXX model with a magnetic field, 
Eq. (|4.2|l . We can now use the BA to compute the 



eigenstates of the Hamiltonian in the critical interval: 
A e [0,2) For values of |A| greater or equal to 2, 

every spin in the system is pointing in the magnetic field 
direction and the state is a product of S'f eigenstates. 
Entanglement has disappeared. 

As implicit in the BA, the ground state has a well 
defined magnetization, S^. This quantity changes ev- 
ery time the tuning of the magnetic field implies a new 
level-crossing. The relation between the magnetization, 
S^, and the ratio A between the magnetic field A and 
the coupling constant J is well defined and can be easily 
obtained from the spectrum of the model. 
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Figure 11: Entropy Sl in a ring of A'^ = 18 spins, 
for blocks of size up to L = N/2 = 9 spins. The 
curbs correspond to the values of the magnetic field, A € 
{0, 0.24, 0.68, 1.05, 1.35, 1.59, 1.77, 1.89, 1.97}, for which there 
is level crossing. For a given L, the entropy remains constant 
in the interval between level crossing, but every time the cou- 
pling is at one of the latter points the entropy value changes. 
The maximum in the entropy is obtained for the antiferro- 
magnetic system without magnetic field, while it goes to zero 
as A approaches 2. 



The Heisenberg model has two limiting behaviors. On 
the one hand, for A = the ground state is antifcrromag- 
netic, with a null angular momentum eigenvalue. On the 
other hand, for A > 2, the ground state of the system 
corresponds to the ferromagnetic state \F). Figs. 111|) 
and H12|l illustrate how the entropy of a finite spin ring 
changes when A varies in the interval [0,2]. Fig. ^ 
shows the way the entropy decreases as A increases, ob- 
taining the zero value when A ^ 2 where the ground 
state turns into a product state. The plotted values of 
A are those for which the ground state has a level cross- 
ing for a ring with iV = 18 sites. The ground state, and 
therefore also its entropy, remain constant in the interval 
between level crossings. Fig. (|12|) plots the way the en- 
tropy increases logarithmically when A decreases, keeping 
the number of traced spins fixed. 
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Figure 12: Entropy Sl in a ring of A'^ = 18 spins, for dif- 
ferent values of the magnetic field A. Stars correspond to the 
reduced density matrix of one spin, Sl^i, and diamonds to 
the reduced density matrix for half ring, 5'i^jv/2=9- 



2. Anisotropic model 



Figure 13; Entropy Sl in a ring of A'' = 18 spins and for 
different values of the anisotropy A. The critical interval is for 
A G [—1, 1], and leads to the upper, superposed curbs. The 
points out of the critical interval are diamonds with A = 1.5, 
stars with A — 2.0 and triangles with A = 2.5. Thus, as we 
change the anisotropy A away from the critical interval, the 
entropy gets saturated by a value that decreases with A. 



Let us now analyze the behavior of entanglement as a 
function of the anisotropy A in a Hcisenbcrg chain with- 
out magnetic field, that is Hxxz with A = in Eq. H4.1|l . 

In this case, the critical interval corresponds to 
A e [—1,1]. For A = 1 we recover the isotropic anti- 
ferromagnetic model, Hxxx- In the case A = — 1 the 
Hamiltonian can be transformed into the isotropic ferro- 
magnetic Hamiltonian, 

i? = -^a/-t?;+i, (4.18) 

by further rotating each second spin by 180 degrees in 
the z direction. Notice that entanglement remains in- 
variant under this transformation. This argument shows 
that entanglement is a robust magnitude and depends on 
fewer details of the system than e.g. its magnetization. 

The anisotropy is a marginal deformation in the in- 
terval A £ [—1,1], in that for any such value of A the 
system has the same large scale behavior. Instead, for 
other values of A a gap appears in the spectrum, giving 
raise to a new length scale. A finite correlation length 
in the system takes over and all correlations decay expo- 
nentially. Fig. H13() shows that in this case the entropy 
Sl gets saturated as a function of L. 



V. CRITICAL VERSUS NON-CRITICAL 
ENTANGLEMENT 

We now turn to analyze the results of the calculations 
described in the previous two sections. 

The big picture emerging from these calculations is 
that there is a clear distinction between the entangle- 
ment in a non-critical chain and that in a critical chain, 
as measured by the entropy Sl oi the reduced density 
matrix pL for L contiguous spins, see Figs. O and l^. 
For all non-critical models, 5*^ reaches a saturation value 



S* as L increases, whereas for critical chains Sl grows 
unboundedly with L. 



A. Non-critical spin chains 

Recall that in the non-critical regime, a spin chain is 
characterized by a gap between the energy of the ground 
state and that of the first excited state. Relatedly, cor- 
relations between increasingly distant spins decay expo- 
nentially, 

(•^f <l) - K) K+l) - cxp(-L/e), (5.1) 

where ^ is the correlation length (we take the distance be- 
tween nearest spins as unit of distance). For non-critical 
chains, we find that the entropy Sl, a growing function 
with the block size i, is upper-bounded by a saturation 
value S* . This value depends on the parameters specify- 
ing the spin chain, and becomes larger as the chain gets 
closer to a critical point or critical phase. For any fixed 
value of the parameters specifying the spin model, Sl 
approaches the saturation value S* for block sizes L of 
the order of the correlation length ^. We can conclude, 
therefore, that the entanglement between a block and the 
rest of an infinite spin chain has a fixed value for blocks 
larger than the correlation length ^. 

Recall from Eq. (|2.15|l that Sl could in principle grow 
as much as L. This means that, for large block size L, 
the entanglement of a block is negligible when compared 
to its maximal possible value, 

S^ 

Van — =0. (5.2) 

L — ^oo L 

Thus, the ground state \^ g) entangles all the spins in 
the chain, either directly or through intermediate spins, 
but the system does only contain a very small amount of 
entanglement at any scale. 
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The above calculations can be used not only to obtain 
Sl, but also to determine the whole spectrum Sp(pl) 
of decreasingly ordered eigenvalues {p\} of the reduced 
density matrix pL for the block (see section IVl C|) . For 
any L, this spectrum contains only a very small number 
of relevant eigenvalues, together with many small, rapidly 
decaying eigenvalues with an insignificant overall weight. 
Thus, in spite of a exponential growth of the rank of pL 
with L, the effective rank (recall section lTl C|) also gets 
saturated, with the saturation value X'^ being reached 
when L is of the order ^. 

From the point of view of numerical calculations, the 
saturation of can be used to interpret the extraor- 
dinary success of the DMRG in non-critical spin chains 
and other one-dimensional systems. Only a fixed, small 
number of eigenvectors of pL must be retained in order 
to capture the relevant degrees of freedom as L increases 
[2^. Similarly, a bounded effective rank X'y as a func- 
tion of L indicates that non-critical spin chains can be 
efficiently simulated using the techniques introduced in 

s. 

B. Critical chains 

The energy spectrum of a critical spin chain is gapless 
and the two-spin correlation function is characterized by 
a power scaling law, 

K<L)-('^r)«L)~i"^ (5-3) 

where q > (there are cases where marginal deforma- 
tions lead to logarithmic corrections to the power law). 
This scaling law implies an infinite correlation length ^. 
For critical spin chains we have obtained that Sl grows 
unboundedly as a function of L. In particular, in all the 
cases we have analyzed the growth of Sl is asymptoti- 
cally given by 

SL^klog^L, (5.4) 

where k is either 1/3 or 1/6 depending on the parameters 
of the critical chain. Thus, for critical spin chains the 
entanglement between a block of spins and the rest of a 
chain grows unboundedly with the size L of the block, 
in sharp contrast with the non-critical case. In the next 
section we will elaborate on the particular shape 1)5. 4|l of 
Sl for critical chains. 

Given fixed computer memory and execution time, the 
DMRG works with significantly smaller accuracy for crit- 
ical spin chains than for non-critical ones. This result can 
now be qualitatively interpreted in terms of the arbitrar- 
ily large entanglement that links a spin block to the rest 
of the chain as L grows. Indeed, a larger value of 5^ 
implies that more eigenstates of pL must be retained in 
order for the DMRG to achieve a similar accuracy as 
when Sl is small (recall Eq. (|2.21|) ). 

We notice that even if Sl for critical chains grows un- 
boundedly with L, we again have that its value is much 



smaller than the upper bound H2.15|l . Sl/L becoming 
negligible for large L. Therefore, even if entanglement 
is the responsible for the appearance of long-range cor- 
relations at a critical point, we find that the amount of 
entanglement in the ground state of a spin chain is sur- 
prisingly small compared to its maximal possible value. 

Finally, in section IVII it will be argued that in d > 2 
dimensions, the entanglement between a d-dimensional 
block of L'^ spins and the rest of a critical spin lattice 
grows as L'^~^. The same scaling is expected to apply 
to non-critical spin lattices. This implies that Xl grows 
exponentially with the size of the block, explaining the 
systematic breakdown of the DMRG in quantum spin 
systems in more than one dimension. 



VI. SCALING OF ENTANGLEMENT AND 
CONFORMAL FIELD THEORY 

The set of results found in the previous sections ex- 
emplify the connection between entanglement and quan- 
tum field theory concepts. The bridge between quan- 
tum information and quantum field theory can be fur- 
ther explored and exploited in both directions. In this 
sense, we shall translate results from black hole entropy 
and effective actions on gravitational backgrounds to the 
language of quantum information. Working in the op- 
posite direction, quantum information natural measures 
of order based on majorization theory seem to find their 
way into renormalization group irreversibility. The cross- 
fertilization of these ideas requires a change of language 
and quite a lot of background on quantum field theory 
that we cannot cover here in a self-contained or satisfac- 
tory way. Nevertheless, we will very briefly review some 
concepts in conformal field theory and c-theorem neces- 
sary to state their spinoff in quantum information theory. 



A. Entropy and conformal field theory 

The study of entropy in systems with an infinite num- 
ber of degrees of freedom has received quite some atten- 
tion in the context of quantum field theory and black hole 
physics. Historically, the thermodynamics of a black hole 
lead to the Bekenstein- Hawking entropy (see for instance 
[s^ l . This entropy is associated to the counting of micro- 
scopic degrees of freedom inside the horizon and scales 
as the area of the black hole. Some authors suggested 
that the origin of this entropy might be rooted in the 
loss of information forced upon an external observer by 
the existence of a horizon. Although this point of view 
is no longer pursued, its analysis comes naturally due 
to the combination of quantum mechanics and general 
relativity which underlies the holographic principle. The 
field theoretical definition of entropy faces the traditional 
problem of renormalization in quantum field theory and 
receives the name of fine-grained entropy as well as geo- 
metric entropy. 



18 



Three computations of genuine quantum field theory 
entropy are worth recalUng. First, Srednicki p5l | consid- 
ered a properly regularized massless bosonic field theory 
in a universe which is divided by an imaginary sphere of 
radius R into its inside and outside parts. He then nu- 
merically constructed the density matrix of the reduced 
outside system and found that its entropy scaled with an 
area law in 3-1-1 dimensions. For l-fl dimensions, the 
scaling behavior of the entropy was found to be logarith- 
mic. Second, Callan and Wilczek put forward the 
concept of geometric entropy in 1-1-1 dimensions. There, 
the power of conformal symmetry was used at full steam 
in order to compute the entropy of a conformal field the- 
ory when reduced to a finite geometry. We shall come 
back to this result shortly. Finally, the third relevant 
computation was carried out by Fiola, Preskill, Stro- 
minger and Trivedi |27j who mapped a regularized field 
theory in 1-1-1 dimension to the Rindler coordinates and 
recovered the logarithmic behavior of the microscopic en- 
tropy. All these computations needed an explicit ultra- 
violet regulator since quantum field theory contains in- 
finitely many degrees of freedom. Note that in our case 
spin chains come equipped with the intrinsic ultraviolet 
cutoff of the lattice spacing. 

Although more restrictive in scope, the computation of 
the geometric entropy of a l-fl dimensional field theory 
brings the advantage that the result is casted in terms 
of the parameters that classify conformal field theories. 
More concretely, the result found by Holzhey, Larsen and 
Wilczek m reads 

SL = '-^\og,L (6.1) 
6 

where c and c are the so called central charges for the 
holomorphic and antiholomorphic sectors of the confor- 
mal field theory. 

Let us briefly recall that conformal field theories 
arc classified by the representations of the conformal 
group in 1-1-1 dimensions. The operators of the the- 
ory fall into a structure of highest weight operators and 
its descendants. Each highest weight operator carries 
some specific scaling dimensions which dictates those of 
its descendants. The operators close an algebra imple- 
mented into the operator product expansion. One op- 
erator is particularly important: the energy-momentum 
tensor T^^. It is convenient to introduce holomorphic 
and antiholomorphic indices defined by the combinations 
T = Tzz and T = Tgg where z ~ +ix^ and z ~ x^^ —ix^ . 
Denoting by |0) the vacuum state, the central charge of 
a conformal field theory is associated to the coefficient of 
the correlator 

(0|T(z)T(0)|0) = ^ (6.2) 

and the analogous result for c in terms of the correlator 
(0|T(z)r(0)|0). A conformal field theory is characterized 
by its central charge, the scaling dimensions and the coef- 
ficients of the operator product expansion. Furthermore, 



unitary theories with c < 1 only exist for discrete values 
of c and are called minimal models. The lowest lying the- 
ory corresponds to c = ^ and represents the universality 
class of a free fermion. 

The central charge plays many roles in a conformal 
field theory. It was introduced above as the coefhcicnt of 
a correlator of energy-momentum tensors, which means 
that it is an observable. The central charge also charac- 
terizes the response of a theory to a modification of the 
background metric where it is defined. Specifically, the 
scale anomaly associated to the lack of scale invariance 
produced by a non-trivial background metric is 

(o|t;|o) = -^i? (6.3) 

where R is the curvature of the background metric. This 
anomaly can also be seen as the emergence of a non- 
local effective action when the field theory modes are 
integrated in a curved background. 

The results we have found in our analysis match per- 
fectly the geometric entropy computation. In the case of 
the XX and Heiscnberg critical spin chains the central 
charge is c = c = 1 and the model falls into the free 
boson universality class. The critical Ising model corre- 
sponds to a free fermion, thus c = c = i. The central 
charge of the theory is seen to play the role of a mea- 
sure of entanglement. The vacuum of a theory of free 
bosons is more entangled than the one corresponding to 
a theory of free fcrmions. Scaling of entanglement is just 
another manifestation of the ubiquitous organizing prin- 
ciple orchestrated by conformal symmetry. The amount 
of surprise (or entropy) obtained when in a given theory a 
new degree of freedom (a new spin in the block) is added 
must follow scaling as dictated by the representation of 
conformal symmetry corresponding to that theory. Be- 
cause the entropy of the reduced density matrix of the 
ground state is not attached to any particular operator, 
it is natural that the central charge is the parameter in 
control of this measure of entanglement. 

Due to the relation between entanglement and the cen- 
tral charge a number of further connections appear. The 
central charge quantifies quantum correlations as well as 
e.g. the trace anomaly. A plausible interpretation can 
now be given in terms of entanglement. The more quan- 
tum correlated the vacuum is, the stronger the break- 
ing of conformal symmetry appears when a curved back- 
ground is present. More relevant is how the theorem on 
irreversibility of renormalization group flows translates 
to entanglement in a way that will be discussed shortly. 

B. Entanglement in higher dimensional systems 

Entanglement in spin chains obeys logarithmic scaling. 
This result was seen to emerge from conformal symmetry. 
It is then possible to apply similar arguments in higher 
dimensions and complete them with standard arguments 
based on the Schmidt decomposition. 
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Let us consider a d + 1 field dimensional theory at its 
critical point. Its ground state is a pure state. Follow- 
ing Srednicki one could now consider the reduced 
system on e.g. a hypersphere 5"*"^ of radius R. The 
division into the interior A and the exterior B of this 
imaginary hypersphere can be used to write the ground 
state using the Schmidt decomposition in terms of pure 
states l^i), associated to each region 



B 



(6.4) 



where pi are positive numbers and the sum ranges up to 
the minimum of the dimensions of the Hilbert spaces for 
A and B. The standard argument follows that both re- 
duced density matrices share the same eigenvalues and, 
thus, the same entropy. Yet, both systems only share the 
hypersurface separating them, so that the leading term 
of the entropy as any infrared cutoff is sent to infinity 
and the ultraviolet cutoff Xuv is sent to zero must scale 
as its "area" 



Cl 



L 



d-1 



(6.5) 



where ci is a known coefficient related to anomalies 
that we shall discuss later. This leading scaling law for 
Xuv also holds for massive particles as checked by 
explicit calculations in free massive theories [s^ . The 
"area" law associated with the hypersurface is under- 
stood as an effect coming from the loss of coherence be- 
tween the points at each side of the boundary separating 
the interior and exterior parts of the universe. It is then 
natural to expect that microscopic condensed matter sys- 
tems will follow the same law. 

At variance with spin chains, the ultraviolet cutoff Xuv 
gets now mixed with the global coefficient ci and it is 
unclear how to extract observable information from the 
latter. It has been shown that ci corresponds to the co- 
efficient of the linear term in the curvature in the effec- 
tive action of a field theory in a non-trivial gravitational 
background |^. Then, ci equals 1/6 for scalar particles 
and 1/12 2^^''^^ for Dirac fermions. More precisely, every 
fermionic component contributes to ci as half a boson. In 
1-|-1 dimensions, we worked with spinless fermions, thus 
the relative factor of 2 between 1/6 for the critical Ising 
model and 1/3 for the XX and Heisenberg critical chains. 
In dimension 3-1-1, the entropy of a system of a free Dirac 
fermion will carry twice more entropy than a free boson 
since the Dirac fermion is made of four components. 

Entanglement is thus also connected through Eq. (|6.5|) 
to the effective action of quantum field theories on grav- 
itational backgrounds. It is a remarkable fact that in 
1 -|- 1 dimensions the effective action has a unique non- 
local form proportional to the central charge c. When 
this non-local action is expanded in powers of the curva- 
ture, all terms carry the dependence in c. It follows that 
the trace anomaly, which is a derivative of this effective 
action with respect to the metric, is also proportional to 



the central charge. This is no longer true in higher di- 
mensions. The effective action of a quantum field theory 
defined on a gravitational background develops infinitely 
many apparently unrelated structures. The entropy of 
entanglement seems to be related to ci 0|, the coef- 
ficient of the linear term in the curvature R (called oi 
in |53]). On the other hand, other contributions to the 
trace anomaly giving rise to the Euler density and to 
non-trivial two-point energy-momentum tensor correla- 
tors are related to the structures that come quadratic 
with the curvature. All the coefficients in the effective 
action seem to quantify different aspects of entanglement. 

A remarkable result concerning the saturation of en- 
tanglement in non-critical systems can also be translated 
from quantum field theory to spin chains. It has been 
proven Q that the leading order result when the ultra- 
violet cutoff is sent to infinity, Eq. H6.5|l is not modified 
by adding a mass to the scalar field. The leading behav- 
ior of entanglement is not affected by going away from 
conformal symmetry because the contribution to the en- 
tropy comes from the entanglement between points near 
each side of the boundary (see also [s^). The scaling law 
in terms of the hypersurface separating them is respected 
although a finite correlation length is present. Moreover, 
the subleading corrections to the "area" law are known 
when Xuv ^ ^ ^ir, Xir being an infrared cutoff which 
defines a separation of space in two regions in analogy to 
our previous R, 
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where the subindex m labels the massive theory. The 
case d = 1 + 1 corresponds to the universality class we 
have been discussin g in systems of spin chains and yields 
the following result [5q|^ 



S" = - - Inmxuv 
6 



1 



m 



(6.7) 



This latter result is indeed observed in our computation 
(and will be presented more extensively in some other 
publication). For instance, the external magnetic field 
deformation of the Ising model (7 = 1,A 7^ 1), Fig. jnj, 
corresponds to an effective mass m ~ |1 — A|. 

The lack of saturation of entanglement in massive theo- 
ries have clear implications for the application of DMRG 
techniques as well as other modifications of the Wilso- 
nian block-spin idea for condensed matter systems. En- 
tanglement always diverges (even off criticality) in more 
than 1-1-1 dimensions and limits the success of such an 
approach. It seems conceptually preferable to construct 



^ We quote this equation from D. Kabat and M. J. Strassler works 
in spite of the fact that the logarithm is in the natural basis. 
Different definitions of entropy are given in different communities 
but the basis just defines the unit of measure. 
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projections of the exact renormalization group based in 
momentum space as the local potential approximation. 
This method has proven quite powerful in scalar theories 
|59l | where it provides good approximations to critical ex- 
ponents in any number of dimensions and can even detect 
triviality. 



C. Majorization and entanglement 
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Let us return to quantum spin chains where entangle- 
ment pervades the system at conformal points. The ab- 
sence of a mass scale makes quantum correlations extend 
to long distances. Consequently, the vacuum structure 
must describe this fact and the entropy shows scaling as 
discussed above. It is arguable that the entropy of the 
vacuum in the reduced system is just one out of many 
possible measures of entanglement. Then, other yet un- 
explored measures of entanglement may bring further in- 
formation about the structure of the vacuum. 

In order to investigate this issue, we can further exploit 
our computations due to the fact that we have explicit 
results for the density matrix of the ground state in the 
Ising and XX models. The 2^ eigenvalues of pz, corre- 
spond to the product of the 2 eigenvalues (1 ± Vm)/^ of 
L fermionic modes, Eq. 



L-l 

n 

m=0 



x„ = 0,1 Vm. (6.8) 



It is convenient to first visualize the typical shape of the 
eigenvalues (l±fm)/2 of the fermionic modes and observe 
that most of the Pl eigenvalues are almost zero or one and 
only a small set will take intermediate values, bringing 
the main contribution to the entropy. Each mode tends 
to remember its pure state origin. 
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Figure 15: Plot of the eigenvalues of the fermionic 

modes, Eq. I3.6UL for the XX model as a function of the 
magnetic field A, Eq. 13. 5L taking L — 200 spins from an in- 
finite ring. Triangles display the eigenvalues for A = 0.5, stars 
for A = 0.999 and diamonds the eigenvalues for A = 0.99999. 
When A approaches 1, the ground state becomes a product 
state and the modes take only the value one or zero. The 
three plots should overlap each other, but we have artificially 
shifted each plot 40 points respect to the previous one. 



The eigenvalues of form a probability distribution. 
It is numerically easy verified that at conformal points 



Pl+2 -< PL 



(6.9) 



that is, the probability distributions associated to the 
density matrices obey a majorization relation. Due to 
the spin structure of the problem the majorization takes 
place at two-spin steps. We have verified that a similar 
result about majorization holds for the quantum Ising 
model at its critical point. 
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Figure 16: Plot of the cumulants of the 8 (triangles), 16 
(stars), 32 (diamonds) and 64 (triangles) larger probabilities 
associated to the ground state pz, as a function of L in steps 
of 2 in the XX model without magnetic field. 



Figure 14: Plot of the eigenvalues 



of the fermionic 



modes, Eq. 13. 601 . for the Ising model, Eq. 13.41 . and L = {50 

(triangles), 100 (stars), 150 (diamonds), 200 (triangles)}. As 
L increases new modes are populated. 



The increase of amount of surprise quantified by the 
entropy is rooted in a deeper sense of ordering in the vac- 
uum. The ground state pz, becomes more and more dis- 
ordered as dictated by a majorization arrow. Note that 
this implies and increasing number of relations between 
the eigenvalues of p^ as L goes to infinity. The vacuum 
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in a field theory is far more ordered than what the scal- 
ing of entropy hints at, in the sense that majorization 
provides a more strict definition of order than entropy. 

In order to establish the role of majorization and its re- 
lation to conformal symmetry, further work is necessary. 
In particular, marginal deformations need to be analyzed 
carefully. Later on we shall argue that a different type of 
majorization seems to hold along renormalization group 
flows. In that case, L is kept fixed but the parameters 
of the hamiltonian change as dictated by renormalization 
group transformations. 



D. Irreversibility of entanglement and the 
c-theorem 

Given the relation between entanglement and confor- 
mal symmetry, some powerful quantum field theory re- 
sults can be borrowed and translated to quantum in- 
formation parlance. It is of particular relevance the so 
called Zamolodchikov's c-theorem that establishes 
that the central charge of a unitary 1-|-1 dimensional 
theory always decreases along renormalization group tra- 
jectories. The existence of a c-theorem in more than 
1-1-1 dimensions has been a subject of a lot of effort 
E0> Ifill- Ifi^ 1^ and a proof has been proposed in 
|65| . Zamolodchikov's theorem in 1+1 dimensions can 
be proven using the spectral densities following Ref. [bJ ■ 
Consider the spectral representation of the correlator for 
two energy-momentum tensors (the energy-momentum 
operator is defined for any quantum field theory and cor- 
responds to a descendant of the identity) 

(o|T„^(x)r^.(o)|o) = 

= I y" c(A, pt) {papp - p^gcf}) {p^Pu - P^9t,u) G{x, A) 

(6.10) 

where A is the spectral parameter (with dimensions of 
mass), c(A, /i) is the spectral function, which depends 
on A and on the subtraction point fi, and G{x, A) is the 
free scalar propagator of a particle with mass A. At a 
fixed point of the renormalization group flow, the spec- 
tral function reduces to a delta distribution c(A, M)lc/t ~ 
c <5(A), where the coefficient c is a constant, reflecting 
the fact that all physical intermediate states are mass- 
less. The UV fixed point can be analyzed taking x ^ 
and it follows that 

cuv = J dXc{X,fi). (6.11) 

where cuv is the central charge of the ultraviolet theory. 
On the other hand, in the IR limit only massless modes 
survive, so that the spectral function can in general be 
written as 

c(A, ^) = ciR (5(A) + Csmootii(A, n) , (6.12) 



where the contribution of all massive modes is contained 
in Csmooth(A, /i). It thus follows that 

Cuv = ciR + J dX Cs„iooth(A, fi), (6.13) 

where the second term on the r.h.s. is necessarily fj,- 
independent. Finally, unitarity guarantees that Cgmooth 
is positive, so 

Cuv > ciR. (6.14) 

This result can be understood as a net decrease of de- 
grees of freedom as weighted by the central charge along 
renormalization group trajectories. Note also these ideas 
match bosonization: in 1-1-1 dimensions two Majorana 
fermions can be made into a boson. All proposed gen- 
eralizations of the central charge in more than 1+1 di- 
mensions arc based on the trace anomaly and share the 
property that fermions weight more than bosons. Exact 
bosonization is then not possible. 

The c-theorem in 1+1 dimensions immediately implies 
that entanglement decreases along renormalization group 
flows. Critical regions in any theory are given by the fixed 
points of the renormalization group transformations and, 
in 1+1 dimensions, they are characterized by the central 
charge of the model. As we have seen, the entangle- 
ment between a block of spins and the rest of the chain 
in a critical region is given by the central charge times 
the logarithm of the length of the block; and from the 
c-theorem, the renormalization flow of any unitary the- 
ory goes from an ultraviolet fixed point to an infrared 
fixed point, i.e., from a higher to a smaller central charge 
value. Then, the entropy of the reduced density matrix 
decreases along renormalization group trajectories, 

Sr > Sl^. (6.15) 

The decoupling of the massive modes in the spectrum 
at long distances irreversibly reduces the amount of non- 
local quantum correlations in the system. This result is 
only valid for unitary theories and is not obvious because 
a renormalization group transformation in quantum field 
theory is made of two steps: integration of high-energy 
modes followed by rescaling. While the first step seems 
to imply irreversibility, the second makes it obscure. Fur- 
thermore, the proper construction of this result implies 
a careful treatment of renormalization. The role of uni- 
tarity is dominant since counterexamples can be build 
where renormalization group flows form limit cycles for 
non-unitary theories. Entanglement loss and unitarity 
are, thus, related. 

There is an explicit result we have already presented 
that can be used to illustrate the loss of entanglement 
along renormalization group trajectories. Consider the 
non-critical massive boson. The departure from the 
massless case corresponds to adding a relevant operator. 
The flow will make the mass grow. The renormalization 
group trajectory is driven by the increase of the mass. It 
must follow from the (l+l)-dimensional c-thcorcm that 
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for increasing mass, the entropy of the system decreases. 
This is indeed the case as verified in Eq. H6.7|l 

5'L,jrii - S'L.ma = -7 In —< mi>TO2. (6.16) 
1712 

A perturbative c-theorem is also known to hold in any 
dimension and, again, it can be checked on Eq. (|6.6(l 



SL.m ~ SL,m=0 < . 



(6.17) 



Furthermore the decrease of the entropy is monotonous 
in m. 

It is natural to try to go one step further and try to 
relate irreversibility of the renormalization group to the 
majorization properties that may structure the entangle- 
ment in the vacuum. Our preliminary numerical results 
seem to indicate that 



PT=° -< PL 



(6.18) 



in spin chains, in all the analytical results shown in 
Eas. H6.6() and (|6.7() and in the computation produced by 
Srednicki |25|]. Although all the above cases correspond 
to renormalization group flows driven by a simple mas- 
sive deformation, one then may speculate that 



(6.19) 



and that majorization is indeed underlying the irre- 
versibility of renormalization group flows. This fact re- 
mains a conjecture, nevertheless there are some numeri- 
cal evidences that will be shown in a later publication. 



in the limiting case of an infinite chain, N — > oo. The 
spectrum of this model was originally computed by Kat- 
sura |23K generalizing the results of Lieb, Schultz and 
Mattis |2l| for the XY model without magnetic field. 



a. Jordan-Wigner Transformation 

The initial spin operators satisfy anticommutation 
rules at any given site but follow commutation rules at 
separate sites. The non-local Jordan-Wigner transfor- 
mation maps these operators into fully anticommuting 
spinless fermions defined by 




(A2) 



{a|, a„i} = Sijn, {ai,a,n} = 0. 



In terms of operators a the above Hamiltonian becomes 

N-l 

1 ^ ^ 



i=- 



(A3) 



+ 7 



ai+iai 



X ^ ajai. 



1 = 



h. Fourier Transformation 
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Appendix A: SPECTRUM OF Hxy 

The material that is shown in this appendix can be 
found in any standard text book that studies spin systems 
(see for instance 0:E3,E3)- Nevertheless, we present it 
to make the article self-contained. In this appendix, we 
determine the spectrum of the XY model with transverse 
magnetic field and open boundary conditions. 



H 



XY 



1 ^ 



2 ^ V 2 

I r, 



-'-^'^ ~x„x , \ In-yn-y 



(Al) 



We can now exploit the (quasi) translational symme- 
try of the system by introducing Fourier transformed 
fermionic operators 



1 ^ 



aie 



(A4) 



-{N - l)/2 < fc < (iV - l)/2. Due to the fact that this 
transformation is unitary, the anticommutation relations 
remain valid 

{dlJpj^Skp V{k,p}. (A5) 
The Hamiltonian now takes an almost diagonal form, 

(N~l)/2 



H 



(_A + cos— )4dfe + 



fc=-(7V-l)/2 
(7V-l)/2 

"2 



(A6) 



27rfc 



~ 1^ sm— (dfed.fc-f 



fe=-(Ar-l)/2 



where an extra term, suppressed by , should be present. 
We shall though ignore it here since our results concern 
the limit N —f oo. 
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iubov Transformation 



A final unitary transformation is needed to cast the 
Hamiltonian into a manifestly free particle theory. This 
so-called Bogoliubov transformation can be expressed as 



bl = Uk dl + ivk d-k 
bk ^ Ukdk~ ivk d}_k, 

where Uk = cos 9k/ 2, Vk = sin 9k /2 for 



COSf^fc = 



-A + cos 



(A-cos^)2+72 sin^ 2^ 



(A7) 



(A8) 



Again, due to unitarity of the Bogoliubov transformation 
the operators {bk} follow the usual anticommutation re- 
lation 

{blbpj^Skp y{k,p}. (A9) 
Finally, the Hamiltonian takes a diagonal form 

(W-l)/2 

H= J2 ^kblbk, (AlO) 

k=-(N-l)/2 

where 



Afc = W(A-cos— j +7'sm2— . (All) 

The thermodynamical limit is obtained by defining = 
2'Kk/N and taking the N ^ oo limit 



H 



271 



A^bib 



with 



Arf, = ~ cos 0)2 -t- 72 gin2 . 



(A12) 



(A13) 



Appendix B: CORRELATION MATRIX FOR THE 
XY MODEL 

In this section we show how to obtain the correlation 
matrix (dma„) of the ground state of the Hamiltonian 



H 



XY 



E 



2 ^ V 2 



1 -^7 



0-210.21+1 



1-7 - . 
— - — 021-1021+2 



-N < m,n < TV - 1. 

These correlators were originally computed by Lieb, 
Schultz and Mattis for the XY model without mag- 
netic field and by Barouch and McCoy (2^ for the XY 
model with magnetic field. 

We start by diagonalizing the above Hamiltonian. This 
can be achieved by means of two canonical transforma- 
tions. Let us define 2A^ auxiliary Majorana operators d 
and e, 



d2k-l 

d2k 

e2k-i 
e2k 



2 

E 



cos ■ 



JV-1 
2 

N-1 
2 

E 

^_ N-l 



2Trkl 
2nkl 



021-1 

021 



021-1 
021 



,(B3) 
,(B4) 



< /c < N/2, that take the Hamiltonian into sum of 
Hamiltonians Hk, 



N/2 



where 



iAk 
4 





HxY = Hk , 








k=0 




d2k-i 


T 


" Cfe 


-Sfe 


e2k-i 




sfe 


Cfe 


d2k 




-Cfe -Sfe 





e2k 




. Sfe -Cfe 






(B5) 



d2k-i 
e2k-i 
d2k 
e2k 
(B6) 

Cfe = cos^fe, Sk = sin 9k- The diagonalization is com- 
pleted by a second transformation that acts indepen- 
dently for each value of k, according to 



b-2k-i 




b-2k 


1 


hk-i 




hk 





Uk 


Vk 


Uk 


-Vk 


Uk 


Vk 


-Uk 


Vk 


Vk 


-Uk 


Vk 


Uk 


-Vk 


Uk 


Vk 


Uk 





d2k-i 




e-2k-i 




d2k 




e2k 



(B7) 

In terms of the 2A^ Majorana operators bp, —N < p < 
N — 1, the Hamiltonian is finally written as 



^^^^J E ^k{b2k-ib2k-b2kb2k-i), (B8) 



as we wanted to show. 

The above two orthogonal transformations define W S 
S0{2N), where 



N-l 



1O2I , 



(Bl) 



-N + 1 <p < N. 



(B9) 



where the Majorana operators d fulfill 

djyi, {CL^yi, Q-n} 25^ 



(B2) 



Then, coefficients gi in Eq. H3.46|l are obtained through 
r-^ = W^r^W, (BIO) 
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where 



e 



1 

-1 



(Bll) 



Appendix C: ANALYTICAL EVALUATION OF 
THE CORRELATION MATRIX FOR THE XY 
MODEL 



b. Ising model in a critical magnetic field.- For the 
Ising model, 7 = 1, and when the magnetic field is in the 
critical point A = 1, the initial expression for gi trans- 
forms to, 



91 



d(t> (—1 + cos (f) cos (j)l — sin sin ( 



TT 



\J {I — COS (pY + sin^ (j) 



-1 



sm ( H — © = — -, TT 



(C7) 



In this appendix we present an analytical expression 
for gi, Eq. H3.47|l . for five particular cases of the XY 
chain. 



d. Ferromagnetic limit 

The ferromagnetic limit corresponds to A 
this case, it is easy to see that 



gi = Su 



V7. 



00. In 
(CI) 



e. Ising model 



For an Ising model, 7 = 1, with magnetic field A we 
have 



1 p"^"^ — A 

For values of A e [—1, +1], 

n _ \pi'l'\-^ — \^ (2m - 1)!! yn i<t>m 



where 

(2m)!! = 2'"m! (2m -1)!! 
obtaining the equation 
m = 



(2m)! 
2™m!' 



(C2) 



(C3) 



(C4) 



(2(m+;+l))! 



(2(m+;))! 



(2'"+i + i(m+/+l)!)2 (2'" + '(m+i)!) = 
(2Tn-l)!! \ 2m+i+l ; ^ n 
(2m)!! ^ ' ^ 

/■r^oo (2(m+i+l))! Y^oo (2(m+;))! 

yl^va=-l-l (2™+' + i(m+/+l)!)2 Z^m=-i (2" + ' (m+i)!)S 

I, (2m)!! — 

(C5) 

For the limits A = or A = 1 these equations reduce 

to 

a. Ising model without magnetic field.- The Ising 
model is recovered for 7 = 1. When no magnetic field 
is present, A = and gi reduces to 



91 = 5- 



(C6) 



/. XX model with, magnetic field 

The XX model corresponds to 7 = with a magnetic 
field in the range given by A £ [—1, 1]. Making a small 
transformation in the general expression, we get 



91 = 



-X + COS0 



TT I A — COS < 



COS l<pd(f>— / cos 







(C8) 



^ sin l(j)^ 1^0 
^ - 1 Z = 0, 



where (jic = arccos(A). In the particular subcase of the 
XX model without magnetic field, A = 0, we obtain 



2 In 
gi = — sm— , 

ITT Z 



(C9) 



which is equivalent to 



gi ~ 0, I £ even 

9. = f (-l)''-)/^ / e odd. ^''''^ 

In 



g. The XY model with critical magnetic field 

For any anisotropy I7I < 1 and critical magnetic field, 
A = 1, the general expression is recast into. 



91 



^ d(j) (cos 0—1) cos l(j) — 'J sin (p sin Icj) 



(1 — cos (/))2 + 72 sin^ • 
sin {I + i)0 



7+1 r 

2n 



^sin^0/2 + 72cos2(/)/2 
7-1 sin(?-i)(/) 



(Cll) 



2n 



sin^ 0/2 + 72 cos2 0/2 



h. The XY model without magnetic field 

Finally, following reference |23| the limit of the XY 
model without magnetic field corresponds to gi expressed 
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as follows: 



where, 



91 =- 
91 =0 
where, 



I G even, 



I e odd 



(C12) 



-,21 I 2' 



1 + 7 



(C15) 



L, 



Tv/2 



\/ cos^ (^ + 72 sin^ (f) 



(C13) 



A series expansion of this integral is given by: 



Lih) = (-1) 
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